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Abstract 

We derive rigorously the leading asymptotics of the so-called Anderson integral 
in the thermodynamic limit for one-dimensional systems. The coefficient of the 
leading term is computed and shown to be strictly positive. This implies a non- 
trivial lower bound for the exponent in Anderson's orthogonality catastrophe 
pertaining to the ground states in systems of non-interacting fermions. 



1. Introduction 



In 1967, P.W. Anderson [2 studied the transition probability between the ground state 
of N free fermions and the ground state of N fermions subject to an exterior potential 
in R 3 . Interestingly, he found that this probability decays like N~~< for some 7 > as 
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1. Introduction 



N — >• oo. Here, we give a rigorous analysis of this so-called orthogonality catastrophe 
for one-dimensional systems. 

To begin with, let us briefly sketch the many-particle problem underlying our consid- 
erations. The state space of N fermions is the TV-fold anti-symmetric tensor product 

H N := H A . . . AH 

N— times 

of some one-particle space % (e.g. % — L 2 {VL) ® (D 2 , il C R d , s,d£N) where a one- 
particle Hamilton operator H : D(H) — > H is defined. Since we assume our particles 
to not interact the corresponding operator H N on H N is simply a sum 

H N ;=H A1A...A1 + -- + 1A...A1AE 

If H has a discrete spectrum consisting of (simple) eigenvalues Ai < A2 < ■ ■ • with cor- 
responding eigenvectors ipi, <f2, ■ ■ ■ one can easily construct the analogous iV-particle 
quantities. In particular, the ground state ip N is a Slater determinant and the eigen- 
value X N a sum, i.e. 

tp N = tpi A . . . A ip N , X N = Ai H h \n- 

Let Hy '■= H + V be a second operator on T~L with eigenvalues fii < /i2 < • • • and 

(simple) eigenvectors ip±, tp2, The operator Hy is defined analogously to H N and 

thus the new ground state is 

ip N = tpi A ... A ^jy, A 1 " = Mi H 1" Mw- 

The transition probability, 2? at, studied by Anderson is given through the scalar prod- 
uct 

V N := \(tp N ,ip N )\ 2 = \det(( i p 3 ,yj k )) hk=1 _ N \ 2 - (1) 
It can be estimated (see Corollary 15. 4p as 

N 00 

V N <e % », In:=Y, E \&jM 2 - (2) 

j=l fe=AT+l 

Here, Tn is the so-called 'Anderson integral' which is the object of our main interest. 
The asymptotics we wish to analyze involves a second parameter L reflecting the 
system length so that H = H L = L 2 ([0, L] d ) is the Hilbert space of (spin-less) fermions 
confined to the box [0,L] d . Therefore, we work with a sequence of Hilbert spaces T-Ll 
and ground states (p N — (p^, ip N — with L > 0. In the thermodynamic limit we 
let N, L — > 00 with the particle density p = N/L d being kept fixed. The main result 
f Theorem 15. 3p is an asymptotic formula for the Anderson integral 

I N L = jlnN + o(lnN), N,L^ 00, 
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2. Representation of the Anderson Integral 



in dimension d = 1 and with a slightly different convention for the box size (namely 
21/ instead of L) and the density p — (N + |)/(2L). We show that the coefficient 7 is 
always positive and strictly positive at least if the potential is small. 

The main ingredient of the proof is an integral formula for In,l ^Proposition 12. 1]) . 
which holds true under rather general conditions. It rests essentially upon the Riesz 
integral formula for spectral projections and Krein's resolvent formula. To adapt it to 
Schrodinger operators we derive a resolvent formula involving abstract differentiation 
and multiplication operators (Proposition 12.21) . Via this formula, a sequence of scalar 
functions comes into play which tends at least informally to a Dirac delta function. 
This is made precise in Sections 13.31 and SI hence the name delta-term and delta- 
estimate. The singularity represented by the delta sequence reflects in a way the 
singular transition from a discrete spectrum to a continuous spectrum as L — > 00. 

Our method requires a rather detailed and precise knowledge of the free Dirichlet 
problem, in particular of the resolvent. Almost everything one needs to know about 
the perturbed problem, however, is can be read off from the so-called T-matrix. The 
perturbed eigenvalues do not enter in the actual asymptotic analysis. We only need to 
make sure that the number of perturbed eigenvalues below some fixed (Fermi) energy 
is asymptotically for large N the same as for the free problem (see Proposition 13.111 
and the discussion before Corollary 15 .4[) . This is related to the spectral shift function 
(see [5] for potentials with compact support). Interestingly, a lot of work has been 
done to derive asymptotic formulae for the perturbed eigenvalues at large energies. 
Except for [T], we are not aware of studies that include also the dependence on L as 
well. 

Anderson's orthogonality catastrophe has attracted a lot of interest in solid state 
physics since its first mentioning. There are early attempts to derive the exact asymp- 
totics of the determinant, 2? at, itself by Rivier and Simanek 11. who used the adiabatic 
theorem. An improved version is due to Hamann [5 who could treat the thermody- 
namic limit only informally. Recent numerical investigations have been carried out by 
Weichselbaum, Miinder, von Delft, Goldstein, and Gefen [TO [5] who also present some 
physical background and refer to further reading. Frank, Lewin, Lieb, and Seiringer 
[4j Eq. (11)] considered the related problem of proving a lower bound to the energy 
difference - in our notation below tr(ifyll — HP) - directly in the thermodynamic 
limit in terms of semi-classical quantities. 



2. Representation of the Anderson Integral 

Let T-L be a Hilbert space with scalar product (•,•), which is anti-linear in the first 
component and linear in the second component and let || • || be the corresponding 
vector norm. The induced operator norm will be denoted with the same symbol || • ||. 
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2. Representation of the Anderson Integral 



We consider a self- adjoint operator H : D(H) H, D(H) C H, and a bounded 
operator V :U^U. Then, iJy = iZ + V is self-adjoint as well with D(H V ) = D{H). 
We denote by cr(H) and cr(Hv) the spectrum of H and i/y, respectively, and by 

R(z) := {zl - H)-\ ze€\ a(H), R v (z) := (zl - H v )-\ z£(C\ a(H v ) (3) 

their resolvents. From spectral theory we know 

= dfatfr.W = dist(,,l(^)) - < 4 > 

We borrow some notation from scattering theory (see e.g. [TH 3.6]). Note, that for 
z (t(H) the operator (1 — VR{z))~ 1 exists and is bounded if and only if z ^ a{Hy)- 
The same holds true for (1 — Hence, the so-called T-matrix 

T(z) := (1 - VR(z))~ 1 V = V(l - R(z)V)- 1 (5) 

exists for z e <D \ (<r(7Z) U <r(i?y)) with 

and is analytic there as a function of z. Krein's resolvent formula 

R v {z) - R{z) = R{z)T{z)R{z) (7) 

relates the resolvents R(z) and Rv(z) with each other whenever the T-matrix exists. 
The operator V plays an important role via its polar decomposition 



V = \V\J = J\V\ = y/\V\Jy/\V\ (8) 

with J* = J = J 2 and || J|| = 1. At times, it is more convenient to work with the 
Birman-Schwinger operator 



K{z) := y/\V\R{zW\V\, z e C \ a(H), (9) 

than with the resolvent. Likewise, instead of the T-matrix the operator 0,(z), defined 
through (cf. m 3.6.1,1]) 



T(z) = y/\V\m(z)y/\V\, Q(z) := (1 - ^\V\R(z) ^V].])- 1 (10) 

is often used. It exists and is bounded for z ^ (D \ (<r(H) U a(Hv))- The converse is 
true, too. That is to say, if z £ c(H) and Q(z) exists and is bounded then z £ a{Hy). 
To see this, one first shows that 1 — R(z)V is injective and has dense range and, in a 
second step, that the range is closed. 
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2. Representation of the Anderson Integral 



2.1. Operators with a common spectral gap 

Riesz's integral formula yields a handy expression for the Anderson integral when the 
operators H and Hy have a common spectral gap. That is to say their spectra can be 
written as 

a{H) = <7!(ff) U a 2 (H), a(H v ) = a^Hy) U a 2 (H v ) (11) 

such that there is a closed contour r C (0 with each o\ being inside and each a 2 outside 
of r. Let P be the spectral projection of H belonging to <Ji(H) and let II be defined 
likewise for Hy ■ The Anderson integral in question is 

X := tr [P(l — II)] . (12) 

In our application P is trace class and hence X < oo. The Riesz formula reads 

P=-L / R(z)dz, n=-!- / R v (z)dz. (13) 
2m J r 2m J T 

Note that both integrals have the same T from above. For our purposes, an infinite 
contour is more appropriate. In particular, due to the special form of the free Green 
function (see (fM|)) a parabola will do best. 

Proposition 2.1. Let P be trace class and assume that the sets o~\^ 2 in (llip are such 
that 

supai(iJ) < v < mia 2 {H) 1 s\vpcji(H v ) < v < mia 2 (H v ) (14) 
with some v € 1R. Then, the Anderson integral ^ 



X = / tr \PR(z)T(z)R(z) 2 T(z)] dz (15) 

2ni '- L J 



with the parabola T u := {z = (^JD + is) 2 \ s 6 R}. 

Proof. First off, P(II - 1) = P(U - P). By Riesz' and Krein's formulae, CE3J) and Q, 

n - P = — / (i? v (z) - R(zj) dz = J- [ R(z)T(z)R(z) dz 
2ni J T 2iri J T 

with the closed contour T from above. Since P is trace class and the other operators 
are bounded we may take the trace. Using the cyclic commutativity we obtain 

X = — [ tr \PR(z)T(z)R(z)] dz = — [ tr \PR(z) 2 T(z)] dz 

2-ki J r 2m Jy 

since P commutes with R(z). Recall that R(z) is differentiable in z € (D \ cr(H) with 
R'(z) = —R(z) 2 . Since all functions involved are analytic for z € © \ (c(H) U a(H\r)) 
we may integrate by parts, 

X = — I tr [PR'(z)T(z)] dz = [ tr [PR(z)T'(z)] dz. 
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2. Representation of the Anderson Integral 



The relation 

T'{z) = T{z)R'{z)T{z) = -T{z)R{zfT{z) 

proves the formula with T. Using the estimates ((H) and © we can bend the closed 
contour into the parabola T u since the integrand decays sufficiently fast at infinity. □ 



2.2. Schrodinger-type operators 

A typical Schrodinger operator is built from differentiation and multiplication opera- 
tors. Let us introduce two operators V and X satisfying 

[V,X] = 1. (16) 

We assume V : D(V) -> H, D(V) C H, to be densely defined on H and X : H -> H 
to be bounded such that XD(V) C -D(V). Thus, ([HI) is meant to hold true on D(V). 
Self-adjointness of Schrodinger operators often results from boundary conditions which 
usually lessen the domain of definition. Let therefore — V 2 have a self- adjoint restric- 
tion H : D(H) H. That is to say D(H) C £>(V 2 ) and 

H = -V 2 on D{H). 

The resolvent of H, 

R(z) = (zl - H)-\ z e C \ <j{H) 
is a well-defined and bounded operator with R(z) : T~L —> D(H). The latter implies 

(zt + V 2 )R(z) = (zl -H)R(z) = 1. (17) 

In general, this equality fails to hold true when the order of terms is switched as can be 
seen in Proposition 13.71 This is the reason why the following resolvent formula gives 
non-trivial results. 

Proposition 2.2. For operators V and X as in (|16t let us assume in addition 
XD(H) C D(H). Then, the decomposition 

R{zf = -R(z) - hxV,R(z)} + -D(z) = ±(R(z) - C(z)) + -D(z) (18) 
z Iz z Iz z 

holds true on _D(V 2 ) and for z € (D \ o~(H). Here, 

D(z) := {^X - R(z)V) [V, R(z)], C(z) := XVR{z) - R{z)VX. (19) 

The operator D(z) is the so-called 'delta-term' and satisfies 

(zl + V 2 )D(z) =0 (20) 

on L>(V 2 ) and for z G <C\a(H). 
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3. One-dimensional Schrddinger Operators 



Proof. We start off from the elementary formula 

R(zf = -R(z) + -R(z)HR(z) (21) 

z z 

and rewrite the last term. By dint of the product rule the commutator in (|18p becomes 

[XV, R(z)] = X[V, R{z)] + [X, R(z)}V. (22) 

Formula JIBJ implies [V 2 , X] = 2V. By noting XD(H) C D(H) and R(z)(zl-H) = 1 
on D(H) we obtain 

[X,i?(z)] = i?(z)[zl - H,X]R(z) = R(z)[V 2 ,X]R(z) = 2R(z)VR(z). 

Thus, 

-[X,R(z)]V = R{z)VR(z)V 

= R{z)V 2 R(z) + R(z)V[R(z), V] 
= -R(z)HR{z) - R(z)V[V,R{z)]. 

We solve for R(z)HR(z), insert this into ([2"Tjl. and use ([22]). Then, 
R( z ) 2 = Ii?( z ) - _L[X, i?(z)]V - ii?(z)V[V, i?(z)] 

= - - hxv, R(z)} + j-X[V, R(z)} - -R(z)V[V, R(z)] 
a Iz Iz z 

= -R(z) - ^-{XV,R(z)] + -D(z), 
z Iz z 



which is the first equality in ([18]) . The second one follows by dint of the commutation 
relation XV = -1 + VX. Finally, by (US]) 

(zl + V 2 )(±X - R(z)V) = ±X(zl + V 2 ). 

Then, 

(zl + V 2 )[V, R(z)} = V(zl + V 2 )R{z) - {zl + V 2 )R(z)V = 
shows ([201). □ 



3. One-dimensional Schrddinger Operators 

We look into the special case of Schrddinger operators with Dirichlet boundary condi- 
tions on the finite interval [— L,L\, Our Hilbert space then becomes R = L 2 [—L,L]. 
Actually, it ought to bear an index L as well as all operators defined on it and related 
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3. One-dimensional Schrodinger Operators 



quantities. However, since this dependence is ubiquitous we tacitly suppress it. In our 
concrete case, 

V = ^, (X<p)(x) = x<p(x). (23) 

The domain D (V) as well as D(V 2 ) can be described with the aid of Sobolev spaces 
which we do not need in detail herein. One can show that X(D(V)) C -D(V). The 
operator H becomes 

H=-V 2 = —^ on D(H), 

where D(H) is Z?(V 2 ) restricted by Dirichlet boundary conditions. Because of that 
we have XD(H) C D(H). The corresponding eigenvalue problem reads 

-ip" = \tp, cp(-L) = = ip(L). 

The eigenvalues Xj and normalized eigenfunctions jjj, j € N, are 



-j= sin(^x) for j even, 
^= cos(^x) for j odd. 



(24) 



We translate the integral formula in Proposition ^ . 1 l and the resolvent formula into 
the framework of Schrodinger operators. For the v G R in Proposition 12.11 separating 
the two parts of the spectrum we choose the so-called Fermi energy 



2L 



(N 



2' 



(25) 



Thereby, the spectrum of H decomposes into 

a(H) = cn(H)Ua 2 (H), ai{H) := {Xj | 1 < j < N}, a 2 (H) := {Xj \ j > N +1} 

and the parabola r„ N becomes what we call Fermi parabola 

r^v := {z — vn — s 2 + H^/vns I — oo < s < oo}, dz = 2i(y / I^v + is) ds. (26) 

The distance of the Fermi parabola from the spectrum is 

\z-Xj\ = \^/v^+is + ^/Tj\\^/v^+is-y/Tj\ > (^7v + s 2 )^((y^- v /A7) 2 + s 2 )^, (27) 

which will be used at various points in particular with s = 0. The spectral projection 
P in the Anderson integral (fT2)l becomes 



N 

The perturbed operator Hy is given by 

H v = H + V, 



(28) 
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3. One-dimensional Schrodinger Operators 



where V is the operator of multiplication by a real- valued function V, the potential, 
denoted by the same symbol for the sake of simplicity. Some results further below 
will be uniform in L. To formulate this conveniently we assume that the potential V 
is already defined on the whole of 1R and not only on the interval [-L, L]. Thus, we 
denote by ||V|| r , 1 < r < oo the i r (R) norms of the function V. Besides integrability 
some statements will require the potential to decay a bit more regularly 

\x 2 W(x)\ < (1+ ^ p )1+a , p > 0, a > 0. (29) 

It is convenient to formulate the special case p = separately 

\V(x)\< (1 + ^ )1+0 . «>0- (30) 

If V € i°°(R) then the operator V is bounded, in concordance with Section [5J regard- 
less of L. In particular, D(Hy) = D{H). Furthermore, since the free eigenfunctions 
are obviously delocalized, V € L 1 (R) implies 

v\ Vj \\ < ^IMIf, (3i) 

which will be used throughout. The spectrum of Hy is given through the corresponding 
Dirichlet problem 

-ip" + Vip = nip, ip(-L) = = ip{L). 

It consists solely of simple eigenvalues, which follows easily via uniqueness results for 
ordinary differential equations. We denote them by k € N with the usual ordering 
Mi < M2 < • • • ■ The decomposition (fTTj) of a(Hy) will be studied in Section The 
normalized eigenfunctions of Hy are ipk, k <= N and the spectral projection II in (|12[) 
reads 

M 

n M :=5^(^,-)Vfc- (32) 



fe=i 



Note that in general M ^ N (see Section 



3.1. Free resolvent 



The spectral representation of the free resolvent ((3]) with ([M]) reads 

oo ^ 

^(*) = Z)— Tfe.>i' (33) 

3=1 ^ 1 

The corresponding kernel or Green function is given by 

_^(sm(y/z(x-L))sm(y/z(y + L)) -L<y<x<L 
[ ' ,V> W{z) | sin(y / i'(x + L)) sm(y/z(y — L)) -L<x<y<L [ ' 
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3. One-dimensional Schrodinger Operators 



with the Wronski determinant 

W(z) = 2v / zsin(Ly / z) cos(Ly/z) = \fz sm(2\/~zL). (35) 

By rewriting the Green function one can cast the resolvent into a form where the L 
dependence is more tangible 

2y/z I sm{Ly/z) cos(Ly>z) * 

The operators P s (z), P c (z), and G(z) have the kernels 

P s (z; x, y) := sin(Vzx) sin( v / zy), P c (z; x, y) := cos(^/zx) cos(Vzy), 
G(z; x, y) := s\n(yfz\x - y\). 

Note that P s (z) and P c (z) are rank-one operators which makes the resolvent differ 
from the operator G(z) by a rank-two perturbation. We would like to apply the delta- 
estimate from Section [4] directly to R(z) and Q(z) (cf. (JTOj) ) . This is not possible 
because of the apparent singularity in the prefactors of P s (z) and P c (z) in (|36[) at 
z = Vff. In a first step we replace z in the benevolent operators P s _ c (z) and G(z) 
by isn and retain the malevolent dependence in the function r. This motivates the 
definition of the operators (i/jy, Ls) and fi* (vn,Ls) in (HT1) and (JHSJ) , respectively. 
In ((54)) '). we estimate the difference between R(z) and R^(vn , Ls). Later, in our 
main Theorem 15.31 we use these operators to compute the coefficient of the leading 
asymptotic ./V-behaviour of the Anderson integral. To begin with, we have a closer 
look at ([311). At the Fermi energy ([25]) 



(39) 



sin(LV^) = ^M) Lf J > cos(L^) = ^C" 1 )^ 1 . ( 38 ) 

which implies on the Fermi parabola (|26p 

sin(L(V^v +is)) = -^(-l) L ^ J (cosh(Ls) + i(-l) N smh(Ls)), 
cos(L(V^v +is)) = -^(-l) L ^ J ((-l) iV cosh(is)-isinh(Ls)). 

Furthermore, we have 

cos,(L(^/vn + is)) J t(Ls) for A even, 
sin(i(y^jv + is)) l-f(La) for A odd, 

where 

. . coshs — isinhs , . , . , . ., m , 

cosn s + % smn s 
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3. One-dimensional Schrodinger Operators 



Now, we keep the s-dependence only in the scalar function r but not in the operators 
P s>c (z) and G(z) and introduce 



1 



2 v /7 7 aF 
1 

2JU^ 



(t{Ls)P s (v n ) - f{Ls)P c (v N ) + G(v N )), 
{-f(Ls)P s (v N ) + t{Ls)P c {v n ) + G{v N )). 



This can be seen, in a way, as the limit of the resolvent as L — > oo (cf. ([5 
that R^^n, s) differs from G(vn) by a rank-two perturbation. 

The operator C(z) in (TT9)) has the kernel (x, y € [— L, L]) 

J xcos(y / i(a; — L)) s\n(y/z(y + L)) + ysm(^/z(x — L)) cos(y / i(y + L)) y < x, 
1 xcos(^/z(x + L)) s\n(^/z(y — L)) + ysm(^/z(x + L)) cos(y/z(y — L)) x < y. 



(41) 



Note, 



(42) 



The Green function and related quantities are to be evaluated on the Fermi parabola 
IV. Most of the crucial estimates needed below come from a bound for W(z). 

Lemma 3.1. There is a constant < Cd < 16 such that for all s S R, L > , N S N , 
and vn as in (|25[) we have 



1 



<C D e- 2L ^, 



| cos(L(y/v^ + is)) \' 



<C D e~ 2L ^. 



(43) 



\sin(L(y/UN + is))\ 2 
Moreover, for z € Tn (see (1261) ) 

\R(z;x,y)\ < C ° e -l«ll»-»l, \C(z;x,y)\ < \c D {\x\ + \y\) e -W*^. (44) 
Proof, (a) For a, s G R, 

| sin(a + is) | 2 = sin 2 a + sinh 2 s. 
We use this in two ways. For \s\ > sq > 0, 



| sin(et + is) | 2 > sinh 2 s > 
For a = ^/vn and all sgR, 



1 + 2 s 



e 2|-l > 



S() 



l + 2s 



„2| a | 



|sin 2 (L(v^ + is))| >sin 2 (Lv^) =sin 2 (-(7V+-)) 



l.s 1 



2" 2 



11 



3. One-dimensional Schrddinger Operators 
We pick an a > to distinguish between large and small s 
\sm(L(^ + is))\ 2 > J [ 



U fork|<f. 
To make the two estimates look alike we use for small s, 

1 = e 2L \ s \ e - 2L \ s \ > e 2L l s l e - 2Q . 

Choosing a — i yields Co — 16 and proves the first part of (I43|) . Starting from 

| cos(a + is) | 2 = cos 2 a + sinh 2 s 
and using cos 2 (L^/i^v) = ^ one can treat the cosine in (|43[) in like manner. 

(b) For L > x > y > -L 

| sin((V^7 + is)(x - L)) sin((V^ + is){y + L))\ < e l«ll*--&l e l«lliH-£| = e \'\<™-\*-v\) } 

where we estimated the sine by the exponential function. For x < y the bound looks 
the same. Using we obtain 

\R(z-x v)\ < - - e l*l(2i-|*-»l) 

' K ' n ~ {v N + s 2 )h 2\sm(L{^ + is))\\cos(L(^ + is))\ 

< Cd . c -\s\\x-y\ 

2{v N + s 2 )^ 

which proves the first estimate in (|44j) . The estimate for C(z;x,y) in ((44)) follows 
likewise. □ 

The optimal a in the proof above is given through the transcendental equation 

a 1 _„ 



1 + 2a y/2 ' 

which can be treated via the Burniston-Siewert method (cf. [3 §14.10]). Solving this 
numerically gives a « 0.8218 and the optimal constant Cd ~ 10.348. 

We study the operators in (|36p and start with the functions uj s {z) 1 u> c {z). 



uj s (z;x) := y/\V(x)\ sm(y/zx), uj c (z,x) := y/\V(x)\ cos{\fzx), igB, Z € <D, (45) 

which come from P s (z) and P c (z) in (j3"T)l . For real values of y^i = a e R, which 
includes the Fermi energy z — vn, one obviously has 

K(a 2 )||<||^||J, || Wc (a 2 )||<||^||f. (46) 
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3. One-dimensional Schrddinger Operators 



The behaviour in the entire complex plane is more complicated. To formulate the 
corresponding estimates in a convenient way we introduce two auxiliary functions (see 
Lemma IA.3[) 

ha( t):=!w^ if *" 1 ' (47) 
lw- ifO<i<l, 

w L , a (s) := -e- L ^ + (2\s\) a h a (Ls) (48) 
a 

with a > 0. Note, that WL, a is continuous at s = 0. 

Lemma 3.2. Let V G L^R) satisfy in addition (|30[) and (|29[) with some a > and 
< p < 1. TTien, 

IK c (z)|| < C> ilQ (*)*e L l*l (49) 
\\u a , e (z) - oj s , c (a 2 )\\ < \a\C u L 1 -"w L , a (s)^e L ^ (50) 

with z — (a + is) 2 , a, s € R, and constants 

C u :=(2C a )i, C' u := (2C a , p )i . 

Proof. To prove (|49]) . we estimate 

\\cj s (z)\\ 2 = [ \V(x)\\sm{{a + is)x)\ 2 dx< [ \ V(x)\e 2 ^ dx. 

J-L J-L 



To prove (|50|) . we hrst note 

| sin((a + is)x) — sin(ax)| = \ix / cos((a + it)x) dt\ < |x||s|e' s:r '. 

Jo 

Hence, 

||cj s (z) - uj s (a 2 )\\ 2 = J \V(x)\ ■ \sin((a + is)x) - sin(ax)\ 2 dx 
< s 2 J L x 2 \V(x)\e 2 ^ dx 

< s 2 L 2(l-p) y L x 2P\ V ( x )\ e ^\ dx . 

In either case we end up with an integral of the form 

L \x 2 W{x)\e 2 \ s *Ux < ( L - f. a * e 2 ^dx < 2C a , p w L , a (s)e 2L ^ , 

J-L V 1 + FlJ 



where we used (|29p and Lemma I A. 3 1 The estimates for w c (z) follow in like manner. □ 
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3. One-dimensional Schrddinger Operators 



The condition on V in (|29|) could be weakened somewhat so as to allow for slightly 
more general potentials. 

Lemma 3.3. Let z = (a + is) 2 , a, s G R. Then, the kernels of the operators P s . c {z) 
and G(z) from (1371) satisfy 

\p s Az;x,y)-P s A^^,y)\<\s\(\x\ + \y\)e^ + ^\ 

\G(z;x,y)-G(a 2 ;x,y)\ < \s\\x - y| e l'H-»l. 

Proof. We write the difference as an integral 
P s (z;x,y) - P s (a 2 ;x,y) 

= / — (sin((a + it)x) sin((a + it)y)j dt 
Jo dt y 

= i I (xcos((a + it)x) sin ((a + it)y) + ys'm((a + it)x) cos((a + it)y)) dt, 



o 

and estimate 

\P,(z;x,y) - P s (a 2 ;x,y)\ < f'\\x\ + \y\)e*^\ + ^ dt < \s\(\x\ + |i/|)e |s|(|x|+|j/|) . 

Jo 

The estimates for P c (z) and G(z) follow in like manner. □ 

One could use Lemmas 13.11 and 13.31 to study R(z) or G(z). However, the applica- 
tions we have in mind require that to be done for the Birman-Schwinger and suchlike 
operators. 

Lemma 3.4. Let V G L 1 (1R) and z G IV. Then, the Birman-Schwinger operator 
K(z) from © satisfies 

\\K(z)\\< ° D \\V\U. (51) 
\Jv N + s £ 

If X 2 V G L X (R) with X as in then for the operator C(z) from we have 

\\V\V\C(z)V\V\\\ < 2C D \\X 2 V\\f\\V\\l (52) 
Furthermore, for the operators G(vn) from (|37[) and i?^(^jv , s) from (|41[) we have 

\\VW\g(vn)V\v\\\ < \\v\u, \\V\v\Ri^N,s)V\v\\\ < ^=nii- (53) 



Finally, if V satisfies (|29[) with < p < 1 and a > then we have 



\W\V\(B*(y N ,L8) - R{{^ + isY)U\V 



< 1 r^=\\V\\ l \s\ + 6Cl P ch 1 - p \s\w L<a (s)e 2L ^, (54) 



where the plus sign is for even N and the minus sign for odd N in i>n . 
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3. One-dimensional Schrddinger Operators 

Proof. Let W be an integral operator with kernel W(x,y) that satisfies 

\W(x,y)\ < \Wi(x)\\f(x,y)\\W 2 (x)\ 

where Wi, W 2 € £ 2 (R) and / e L°°(R 2 ). Let 9? G L 2 [-£, L]. By the Cauchy-Schwarz 
inequality 



— L 



W<p\\ J = \Wi (x) V / /(x,y)W 2 (y)^(y)d7/ 



2 



L 

L r L 

2 / I JV„. .ATI/ ("™.M2 J„, J„.ll,„l|2 



<ll/llLlM2«ilM| 2 , 

where the norms of W\p and / pertain to 1R and R 2 , respectively. Hence 

HWH < H/HoollWillallWalla. 



To prove ([5T1) we can take W\ ; 2 = y\V\ and / = 1 by dint of By the same 

estimate we can prove 1(52]) by using Wi = X y / \V\, W 2 = y/W~\i / = 1 an d Wi = vT^L 

w 2 = jcVJTT, / = 1. 



Because of the obvious bound \G(i>N;x,y)\ < 1 (cf. ((37)) ) we obtain 

\\vw\g(vn)VW\\\ < mk, 



which proves Q for G(v N ). From (gBJ one infers || y/\7\P a ,c{vN) Vrlll < ll^lli- 
Along with |r(Ls)| = 1 this implies for all s S R 



viJ^K.ajviviii 
< ^=(llv^^(^)v^ll + \\VW\Pc{vn)VW\\\ + \\V\v\g(vn)V\v\\\) 



<^=\\vh, 

lJv N 



which is ([53")) for R^iy^.s). 

To prove (|54[) let V satisfy (|29l) with < p < 1 and a > 0. Using the kernel estimates 
in Lemma l3. 31 along with Lemma lA.31 we obtain 

\\V\V\( p :c(*) ~ Ps,cM)y/\V\\\ < 2Cl P ch 1 - p \s\w L , a (s)e 2L ^, 
\\V\V\(G(z) - G(v N ))y/\V\\\ < 2ci P CU 1 - p \s\w L A S )e 2L \ s \ 
which would also be true for other real values than vjq but that is not needed here. 
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3. One-dimensional Schrddinger Operators 



Then for even N 



R+(u N ,Ls)-R(z) 
1 



2(v^v + is) 

1 1 
+ 2 { V^ 



(Ls)(P a {v N ) - P s (z)) - f{Ls){P c {v N ) - P c (z)) + G{v N ) - G{z)] 



1 



-{Ls)P s {v N ) - f{Ls)P c (v N ) + G(v N )] . 



/Vfq + is 

For odd N replace r by — f. In either case we obtain 



□ 



3.2. Truncated free resolvent 



Let Sn{z) '■= PnR(z) be the truncated resolvent with the spectral projection from 
(|28p. We need to control Sn(z) on the entire Fermi parabola (see (p26")) ) and, with 
more care, at the Fermi energy (|25l) . 



Lemma 3.5. Let V G L 1 (1R) and z £ F/v- Then, with Sn(z) = PnR(z) we have 

' ln(iY + l), 



IVISjvWVI^I < -ll^l 



7T + -S 2 

ice r 2 

|v|(s w (*) - SwfaOVNI < —\\v\\i\8\— - 



(55) 



(56) 



Proof. We start off from the spectral representation of Sn, 

N -, 

\V\S N (z)VW\ = Y l 7^(V\V\<Pj,-)VW< Pj . 



3=1 



Applying the estimate ([51]) and then using ([?f| we obtain 



— 1 " 1 2 1*1 



Likewise, using 



z - v N 



4v N + s 2 



< 4- 
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3. One-dimensional Schrddinger Operators 



we find 

JV 



\\VW\(Sn(z) - S N (u N ))V\V\\\ <j\\V\\iJ2 i T n^ 1 x i 

L f-j \z- Xj\\v N - Xj\ 



o 

< 2|s ' 



1*1 1 
IjWVh^-— =- (58) 



3- 



N 



Applying ^ to the sums in (|57l) and (|58[) we obtain (|55[) and (|56p . □ 

The asymptotic analysis in Section r5.2l is based upon a formula for the kernel Sn(z; x, y) 
of the truncated resolvent. In particular it helps to separate the x, y and N dependence 
for special real values of z including the Fermi energy (|25[) . 

Proposition 3.6. Let z > such that i/z > SrN . Then, the kernel SN(z;x,y) of 
the operator Sn(z) = PnR{z) decomposes into 

S N (z; x, y) = S ,n(z; x, y) - Sx,n{z; x, y) - S ,n(z; x, y) + §x >N (z; x, y) 

with 

^ rz,~ ^ f°° „-^^ sinh ((^ + 



S ,n(z; x, y) := 2 ^^- cos(Vz(a; - y)), x N := J e - Vzv — 2 dv, 

k N ^_ , - -i\N f°° „-^^ cosh((iV+ %) v) 



S , N (z;x,y) := ^=cos(Vi(x + y)), A* := (-1)" ^ e -^' u cosh — 2 ^ 



an 



d 



S hN (z; x, y) := ^ J q sm (- Vi(u - -2^-)) ~u **, 

S hN (z;x,y):=— Tz {-1) Bm {—^ u -—— ) ) du. 
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3. One-dimensional Schrodinger Operators 



Proof. With the eigenfunctions from (|24|) we can write 

S N (z;x,y) 



N 1 



l<j<AT 
/ even 

JV 



1 .cosf' 3 * 1 -"'' 



Z — A., 



2L 



N 

E 



l<j<AT 
j odd 

( - 1)J ' C0S r 7rj(a:+y)< 



z — A., 
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where we used the product formulae for sine and cosine. To sum the series we write 



the fraction as a Laplace transform. It is convenient to put z — 



4L 2Z ■ 



Then, 



_ 4L 2 1 , 1 1 \ 4L 2 1 



Z — Xj IT 1z Z — j Z + j 7T 2 Z Jq 

since z > j by assumption. Hence, 

AT 



1 cosh(jv) 



. 2X1 
b N {z;x,y) = -j- 



7T^ 2 







e 2:11 E C0S Mi? 



,irj{x-y). ,%j(x + y)s 

COS ( — j - ( — 1) cos (- 



2L 



2L 



dv 



and with the aid of the addition theorem for the cosine we obtain 

_2 



7i~ 



zS N {z;x,y) 



= 2 



Re 



N 

£ [ cos b'( 



- y) 

2L 



+ iv)] — (— iy cos [j( 



Re 



sin [ (jV +!)(£(£_!) +w) ] 

sin[ i ( zfa) +iu)] 



2L 



_ ( i)Jv cos[(iV+|)(^l + ^)] 



cos[i( 



1 / 7r(x+y) 

2 V 2L 



TO)] 



de 



=: Re/ S - (-l) iv Re/, 



Here we used the summation formulae for the cosine series. The integral I s can be 
evaluated by a change of the integration contour. To this end, let 

ir(x — y) , , „ T 1 
w = a + iv, a := y , M := N+ -. 

We take T = T x U T 2 U T 3 U T 4 , 

I\ = {a + I < v < R}, T 2 = {u + iR | < u < a}, 

r 3 = {»« | o < « < i?}, r 4 = {u I o < u < a}, 
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3. One-dimensional Schrddinger Operators 



orientated counterclockwise. By Cauchy's integral theorem 

^MMwl dw = 
sin -f 

since the integrand has only removable singularities. Because of 

jzw ^jMw) _ oi z {u+iR) sin(M (u + iR)) , g _ M+ i )R „ 
sinf sin(±(u + iR)) 

and z — M + ^ > the integral over vanishes as R — > oo. Hence, 



e lza / s = / e^°+™> sin ( M ( a + iv )) dv 
sin(|(a + iv)) 



and furthermore 



sin^w) 

/- n f" -^sinh(Mu) T . ..sin(Afu) , 

Re /, = cosfza) / e , - dv — / sin z «-a — du. 

Jo smn f Jo sm | 

This gives the terms So.n and Si, at- The integral I c can be treated in like manner and 
the proof is finished. □ 



3.3. One-dimensional delta-term, D(z) 

The delta-term being not trivial, reflects on an abstract level the boundary conditions 
used in the definition of H , which make up the difference between H and — V 2 . 

Proposition 3.7. Let z € <D \ cr(H). For <p € D(V 2 ), and R(z) the resolvent of H 
we have 

(R(z)(zl + V 2 )^)(x) = - Sin( ^ ± L)ML) - p r f {x ~ L)){P{ - L) . (59) 
Furthermore, the delta-term D{z) from (j!9[) reads 



k( I P(z) + L 

4 V s in 2 (yiL) ' y ' cos 2 (^zLY 



D{z) = -{ _. , — P a ( g ) + — ^ P c (z)) (60) 



with the projection operators P s (z) o/nd P c {z) from ((37]). 
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3. One-dimensional Schrddinger Operators 

Proof, (a) To derive ((59)) we integrate by parts two times using the Dirichlet boundary 
conditions, i.e. R(z; x, ±L) = 0, in the first step 

(R(z)V 2 tp)(x) 

i-L 

R(z;x,y)ip"(y)dy 

L 

[R(z; x, yW{y)t L £ dR ^ y \ \ y ) dy 

v{y) dy 



8R(z;x,y) A L f L d 2 R(z;x,y) 



dy 



My) 



-L 



dy 2 



dR(z; x, L) + Wvp-Q _ z f L ^ dy + 



dy ' j_ L 
From the explicit form (|34|) of R(z;x,y) we deduce 

oy sw{2^/zL) 

which implies (|59l) . 

(b) To obtain Formula (|BDJ| we use ^ along with (|5T)|) . 

0= (i?(z)(zl + V 2 ).D(z))(x,y) 

n , n sin(yi(x + L)) . sin(Vi(x - £)) 

= D{z;x,y) - - . " D(z;L,y) + - . D{z;-L,y). 

2 sm(-^/zi) cos^zi) 2 sm(y'zL) cos(-^/zi) 

By the Dirichlet boundary conditions 

(i?(z)V[V, ii(*)])(±L, y) = J R(z; ±L, t/)(V[V, i?(z)])(y', y) dy' = 0. 
Hence, by definition (ITT))) and the explicit form of -R(z) we get 

U(*5 ±L, y) = ±^ x R(z; ±L, y) = ^ ± L)) 



2 y '■ ,! " 2 2 sin(^L) cos(yzi) ' 

Putting everything together we obtain 

L sin(y / i'(a; + L)) sin(y / z(y + L)) + sm(y/z(x — L)) sin(v / i(y — L)) 

D(z\x,y) = = — — — , 

y ' ,yj 8 sin 2 ( A /ii)cos 2 (V5i) 

which implies the statement via the usual trigonometric formulae. □ 
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3. One-dimensional Schrodinger Operators 



3.4. Perturbed resolvent 

Since the perturbed operator enters only through the T-matrix and the operator Q(z) 
(cf. ^ and (flU|) ) we have a closer look at suchlike operators. Recall from Section [5] 
that we already know those operators to exist for z ^ H. What is new herein is the 
uniform bounds on the entire Fermi parabola Tn including vm- 

Lemma 3.8. Let V G L 1 (1R) and assume in addition 

g:=-^=||Vli<l. (61) 

Then, the operator f2(z) from (|10l) exists for all z £ Tjv and is uniformly bounded with 

\Mz)\\ < — ^- =: C n . (62) 
1-9 

In particular, ^ a{Hy)- If in addition V € i°°(lR) then the T-matrix, T(z) (see 
exists and is bounded with ||T(z)|| < Cnll^Hoo. 

Proof. For the operator in ([TUf we have 

\\V\V\R(z)VW\J\\ < \\K(z)\\ < -^j=||V||i < 1. 

For the Birman-Schwinger operator we used (|51[) in Lemma l3.4l Hence, by a Neumann 
series argument we see that fl(z) exists and is bounded with (|62p . Since fjsr ^ °~(H) 
by construction the remark following (|10p shows i/jv £ cr(-ffy). Furthermore, 



HT(z)|| = \w\v\m(zW\v\\\ < WmuWJWWmW = IMUIfiOOII 

completes the proof. □ 

In Section [3TT1 we argued that we prefer to replace the resolvent R(z) by R^Vff , s) . 
Likewise, we replace the operators fi(z) by operators fi* (v^, s), 



fi± (u N , s) := (1 - y/\V\Rt(vN, s)^/\V\J)-\ 

fi<»,oM ■■= (i - ^-^VWigmVWij)- 1 , 

IJV N 



where i?^,(z^v, s) is as in (BT1) . 
Lemma 3.9. Let V £ L X (]R). // 



||V|| 1 <1, <7oo,o := t^HI^IIi < !> ( 64 ) 
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3. One-dimensional Schrodinger Operators 



then the operators f2^(^/v,s) and ilooflivN) (defined in (|63[) ) exist and are bounded 
with 

||^(w)|| < y^— =: C Qoo , ||noo,o(^)|| < — ^ =: C noo , . (65) 

// in addition V satisfies (|29[) a > and < p < 1 i/ien f see (|10p and (|48p _) 

< (7^1^1(1 + ^-"tfli.aWe"!'!), (66) 

where z = (^Jvn + is) 2 and 0+ (z//V) ^^d^N, Ls) is for even and odd N, respec- 
tively. The constant (see ([62]). p9 ]) . (|30 | ) 

C[j :=ma X {-4=||^||i,6CnCn oe< C'l, p cl}. 



Proof. We know from that 



Thus, a Neumann series argument shows that f2^(z//v,s) exists and is bounded with 
the bound given in (1651) . The operator f2 00 ,o( I/ Jv) can be treated in like manner. For 
([61)]) note 

- ft^>iv,£s) - r!(z)v^«(^,L S ) - fl(z))^7|fi± (v N ,s). 
Hence, (1521 . (IM|) . and the first part (|55|) conclude the proof. □ 

We had seen in (|36[) that i?(z) differs by a rank-two perturbation from some operator 
G{z). Thus, one could think of applying the finite dimensional Krein formula to 
fl(z). Unfortunately, the operator 1 — y/\V \G(z)y/\V\ might not be invertible. At 
least, if it should be this cannot be seen by a Neumann series argument because the 
norm estimate would contain an exponential factor as in Lemma 13.21 This is why we 
introduced the operator J7^,(z/jv, s). We only need to know how it acts on the functions 
^s,c{vn)- To this end, we introduce the 2x2 matrices 

o±/„ a \ ._ (Mm), JQ,t>( v N,s)u} s (v N )) (u s (is N ),jn^(v N ,s)u} c (is N ))\ 

°° l N ' >- \(ucM,Jitt(yN,8)u t (v N )) (u c (v N ),JSl±(v N) s)u c (y N ))J ' [b '> 
and 

N) ' \(u} c (v N ),JQ 00 fi(l'N)u s (l' N )) (u} c (i / n),J^oo,o( 1/ n)^c(i'n))J' 

Note that £Ioo,o(vn)* = ^oo,o(^jv) since Jfi o,o( J/ Af) is self-adjoint. 
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3. One-dimensional Schrodinger Operators 

Lemma 3.10. Let V E L 1 (R) satisfy in addition (see (|65|) ) 

1 



iv N 

Then, the 2x2 matrices ^(u^^s), 



ll^lliCn«».o < !• (69) 



Z + (v N ,s) := (1 - -l=Cl+ (u n )t{8))-\ 
IJv n 



Z-{v N ,s) ,= ( 1 + -±—h^ Q {v N )T{s)*)- 1 
lJv N 



(70) 



exist. Here (for t see (|40[) ) 

f(s) := diag(r(s), -f (a)), f (s)*r(s) = 1. (71) 
TTie matrix elements of Q^^visr , s) are given through 

Cl+(m,s) = 2^(Z+(u N ,s) ~ l)t(s)* = Z+is)^^), 



Q^vn^) = -2jv^(Z (v N ,s) - 1)t(s) = Z (^v,s)f2oo,o(^jv), 
where the plus sign is for even N in and the minus sign for odd N. 



(72) 



( Cll 


Cl2 \ 




\C 2 1 


C22/ 





Proof. For operators having the form of fiooC^iVj s) simple algebra shows the formula 

(A - oi(/i, - oa(/ 2j Oss)" 1 = A" 1 + ai [cii(/i, A" 1 -) + da(/ 2 , A" 1 -)] A^ 

+ a 2 [c 2 i(/i,A- 1 -)+c 22 (/ 2 ,A- 1 -)]A- 1 52 

with an invertible operator A, vectors /i )2) 31,2, and ai j2 G C The coefficients are 
given via 

l-a^A-V) -a 2 (/i, A -1 ^) \ _1 
-a^/a.A- 1 ^) l-a 2 (/2,A- 1 32 )y 

To evaluate the inverse on the special vectors $1,2 we recall the definition of the Cjk 
and obtain 

(A - ai(/i, -).9i - a 2 (/ 2 , Os^'V = cnA _1 gi + — c 2 iA~ 1 g 2 , 

ai 

(A - ai(/i, -)gi - a 2 (/ 2 , ■)g 2 )~ 1 g 2 = —ci 2 A~ 1 g 1 c 22 A~ 1 g 2 . 

a 2 

The first matrix element is 

(A, (A - a x (/i, -)3i ~ a 2 {f 2 , -)g2)~ X 9i) = cu(/i, A"V) + — c 2 i(/i, A _1 p 2 ). 

01 
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3. One-dimensional Schrodinger Operators 



Once again, we can compute this by using the definition of Cjk- The remaining matrix 
elements are computed in like manner. Now, for even N put 

7~(s) 7~(s) 

a i = n ; — , «2 = - „ , — , A, 2 = J*u SyC {v N ), g 1>2 = u> SyC {v N ), 
2.Jv N lJv N 



to obtain 

^oo(^jv,s)w s (^at) = cii(s)fi 00)0 (&'jv)w a (i'jv) - f(s) 2 c 2 i(s)il oofi (v N )u> c (v N ), 
tloc(vN,s)uj c (v N ) = -t(s) 2 ci 2 (s)r2oo i0 (v n )lu s (is n ) +c 2 2(s)n (X j,o(i / Jv)w c (i/jv)- 

With Z + (i/jv,s) := (cjfc(s))j,fc=i,2 this gives the matrix elements (IT2"j) . To show that 
is well-defined we look at the entries 

\(u) s>c (u N ), JQoo,o( v n)Us,c(vn))\ < ||^||i||^oo,o(^Jv)||- 
We thus have the estimate 

' M Aoo,o(^)r(«)|| < -4=||V||i||fioo,o(^)|| < -^=||tliCn D0 .o < 1, 



which proves the statement via a Neumann series argument. For odd N replace ev- 
erywhere t by — t. □ 



3.5. Perturbed eigenvalues 

One important consequence of Lemma 13.81 is that the spectrum cr(Hv) of the operator 
Hy on L 2 [— L, L] can be decomposed with respect to the Fermi energy i/jy, 



(73) 



a(H v ) = a 1 (H v )Ua 2 (H v ), 
<ti{H v ) := | fij < i/ N }, a 2 (H v ) := \ ^ > v N }, 

Equivalently, there is an M = M(N) with 

Hj < u N , j = 1, . . . , M, Hj > "N, j > M + 1. (74) 

Exactly N free eigenvalues lie below v^. We need to know how many perturbed 
eigenvalues do so which amounts to estimating M. We will restrict ourselves to non- 
positive potentials. The reason is that for such potentials we have the immediate lower 
bound M > N since the eigenvalues move leftward. For the upper bound we modify 
Bargmann's inequality on negative eigenvalues (cf. jTUJ Thm. XIII.9]). 

Proposition 3.11. Let V <G L°°(1R) be non-positive, V = —\V\, and satisfy (|3H)) with 
a > 0. Let furthermore E > 0. Then, 

or 

M := <E}< —VE + 1 + C E , 

C E :=^ + 1(\\V\\ 00 +E^&^^. 
E 7T " " ' K 2E J a + l 
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3. One-dimensional Schrodinger Operators 



In particular, with E = being the Fermi energy the bound becomes 

M < N +^ + C E - (76) 

Proof. By a shift of the spectrum M equals the number of negative eigenvalues of 

-V" - (\V\ + E)i> = flip, fj(~L) = = 

The eigenfunction ipM corresponding to flu has exactly M + 1 roots. It is convenient 
to distinguish between negative and non-negative roots. We start with the latter 

< xq < x\ < ■ ■ ■ < x m = L 

assuming there are to + 1 of them. Then, fiM is a negative eigenvalue for the Dirichlet 
problem on each Ik := [xk, Xk+i], k — 0, . . . , m — 1. We want a lower bound for the 
distance of two consecutive roots. To this end, we estimate 

(\V{x)\+E)\fj M (x)\ 2 dx<snp(\V{x)\+E) [ \^ M {x)\ 2 dx 

< snp(\V(x)\+E)( Xk+1 ' Xk ) 2 [ \^' M (x)\ 2 dx, 

x£l k K Jl k 

where we used Wirtinger's inequality (see [3]) or in other words, the variational prin- 
ciple for the lowest Dirichlet eigenvalue. If we had 

sup(\V(x)\+E)( Xk+1 - Xk f<l 

then we would conclude with the aid of the differential equation and the Dirichlet 
conditions that 

(\V(x)\+E)\Mx)\ 2 < f W M (x)\ 2 dx 

Ik J Ik 

(\V(x)\ + E)\tP(x)\ 2 dx + fi M I \^m{x)\ 2 dx. 
h -lik 

This is impossible for fiM < and hence, 

1< {Xk+1 ~ Xk)2 sup(\V(x)\+E). (77) 

n x(Llk 

We proceed in two steps. At first, we derive a rough but uniform bound by using 
in (77]). This yields 

Xk+i — Xk > r =: $ an d Xk > kS, < k < m. 

{W\\oO+E)2 
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4. Delta-estimate 



In the second step we use the decay ((30)) of V in (|77[). 

* " " (1 + ^ - (Xfc+1 - + 2(1+S)^ 

This can be cast into the form 



1 < (Xfe+i - Xk) + 



it 2E(l + k5)^ . 1E 



_ ( , -u£ 

Summing up from to to — 1 we obtain 



m< L -\ > tt- < L H ■ 

" 7T 2£ ^ (1 + + % - 7T 2£ <T2.E' 1 + a 

Here we used that the series converges for a > 1 and estimated the sum by an integral. 
The bound for the number to' of negative roots can be derived in like manner and looks 
the same. Adding both inequalities and noting to + 1 + m' + 1 = M + 1 we prove (|75p. 
Finally, ([76]) is an immediate consequence of the definition of v^. □ 



4. Delta-estimate 



Proper treatment of the dominant term hinges on certain integrals of the function 
w a ^L- More precisely, assume we are given a family of functions fi : R + — ¥ R + , 
L > 0, and rate functions O : 1R + — > R + such that 

h(s) < sO(L) and f L {s) < d(L), L > 0, s e R + . (78) 

Proposition 4.1. Let the functions Jl : R — >• R + satisfy ([75)1 wii/i aLd(L) > Q(L) 
where a > 0. Let furthermore a > 0. Then, with WL.a defined in (|48p we /iowe 

Proof. Due to the different behaviour of the two functions making up u>l iQ the proof 
is split up into two parts. 

(a) The exponential function needs only the Lipschitz continuity in (I78[) 

L r° y 1 e - £ '/ L ( a )ds<e(L)- r se - Ls ds^-^. 

Jo Va 2 + s 2 a J a L 



26 



4. Delta-estimate 



(b) Since h a does not decay fast enough we have to proceed more carefully. Pick some 
i] > \ to distinguish between the different regimes in (j78|) . Thus, 



L / 7 9 „ fL(s)h a (Ls)ds 
Jo Va +s 

<-e(L) s a+1 h a (Ls)ds + LfflL) hJLs)ds 
a Jo Jr, Va 2 + s 2 

16(1) /- ir ' Q+1 , , w f , . , 

s a+L h a (s)ds + — — / -====/i a (s)ds 



- a L a 1 a L a Jn. y/l + s 2 s 

<-^[e(L)r? + ^(L)ln(l + -)]. 
a L a 1 rj J 

Recalling the assumption aLi9(L) > O(L) we may choose 

0(L) > 1 



rj = a- 



J &(L) ~ V 

which proves the second part of the estimate. □ 
We also need a similar estimate without the factor f^. 

Lemma 4.2. Let a > and a > 0. Then, with WL, a defined in (|48[) we /iowe 



f 00 1 111 

L / u- L , a ( S )ri S < -{- + — [l + 721n(l + Q L)]}. 

Jo Va 2 + s 2 a a L 

Proof. Once again, we treat the two functions in ici, a separately, 
(a) For the exponential function we have 

' -.e- Ls ds< -L ( e~ Ls ds = -, 



la yja 2 + s 2 a J a 

which proves the first part. 

(b) Once again, the function h a requires more care. We have to distinguish between 
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5. Asymptotics 



small and large arguments. Thus, 

f°° s a 1 fi f°° s a 

L h a (Ls)ds < -L s a h a {Ls)ds + L h a (Ls)ds 

Jo Va 2 + s 2 a J Ji sja 2 + s 2 



11 1 f°° 1 1 , 

:— as 



a L a L a Ji Va 2 + s 2 s 

11, f aL 1 \ 
= --r-i 1 + / ds) 

< i-^(l + V21n(l + aL)). 
a lj 

This concludes the proof. □ 



5. Asymptotics 



In the thermodynamic limit the particle density p is kept constant. Usually, that would 
be wr. However, taking 

N + - 

p:=^>0 (80) 

will make our formulae handier since vn = Tr 2 p 2 is constant then. We start off from 
combining Propositions 12.11 12.21 and 13.71 and write 

tr[Pjv(l-n M )] = -L / \^-tiA N (z)dz + ~ixB N (z)] dz, (81) 

where M = M{N) according to the decomposition ([75)1 . ([71)1 of the spectrum a(Hv)- 
The Fermi parabola Tjv is defined in ([26| . The operators in (|8Tj) are 



A N (z) := P N R{z)T{z)(R(z) - C(z))T(z), B N (z) := P N R(z)T(z)D{z)T(z), (82) 

with the operators Pjy, R(z), C(z) defined in l[2"8"|). (|3"3")l . ([TO)) , respectively. The traces 
can be treated further. For An(z) we use the </jj's and write 



tvA N (z) = J2 —^-(v 3 ,T(z)(R(z) - C(z))T(z) Vj ). (83) 

For Bisf(z) we recall the projection like form of P s (z) and -P c (z), see (pJTl) . Then, 

trB w (z) = d s (2)(J!](^ s (z), ^S N (z)#J!l(z) Ws ( Z )) 

/ / 84 

+ dc(z)(jn(z)u c (z), ^/\V\S N (z)y/\V\JSl(z)u c (z)), 
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where we abbreviated 

ds{z) 2sin 2 (Vi)' dc{Z) ' = 2cos 2 (Lv^) 

with oj s c (z) as in (|45|). The complex conjugates in (|84l) are due to the sesquilinearity 
of the scalar product. We will see that both tr Am{z) and tr Bn(z) decay sufficiently 
fast on the Fermi parabola such that the integrals can be evaluated separately. 



5.1. Subdominant term 



We discuss the subdominant term arising directly from the integral formula. Additional 
corrections will appear in Section 15.21 



Proposition 5.1. Let V € L 1 (R) satisfy (j6Tj) . Furthermore, assume X 2 V € L 1 (M,) 
with the operator X from (|2"5|) . Then, there is a constant C su b > such that 



2z 



tr An(z) dz 



_i ^sub ~ I "ft 



N 



(85) 



where the integral converges absolutely. The operator An is defined in (|83[) and the 
Fermi parabola Tn is defined in (J2HJ) - 



Proof. We start by bounding tiA^(z) for z £ Tjv, i.e. \pz — ^Jv~n + is. From (l83l) we 
obtain 



N 

trA N 



N 

The matrix elements can be estimated with the aid of Lemmas 13.41 and 13.81 as follows 



\{ Vj ,T(z)R(z)T(z) Vj )\ = |(yWI^(2)#(z)WVI^)I 
<||^^ 2 ||Q(z)|| 2 ||i^)|| 
Ci 1 



< 



'-■ V " \ ! 

and 



\{^,T{z)C(z)T(z) Vj )\ < \W\V\vA\ 2 \m\\ 2 \\VW\C(z)V\V\\\ < f 
with constants 

C x := \\V\\ 2 C 2 C D , C 2 := 2||V||f \\X 2 V\\lC 2 n C D . 
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To treat the remaining sum in (|86[) we bound (|27[) from below via 

\z-\j\ > {v N + s 2 )^{{^~ V^f+s 2 )^ > V2(v N + s 2 )^ - v/A7)2v^. 

With the aid of EH we obtain 



N 1 
51 u _ 



< 



1 



1 \/2ly> 



1 



< 



4 ,/L(iV + l) 



for s ^ and thus 

Ci 



1 



4 VTvTT 
A (fjv + s 2 )* VT«T ^("n + s 2 )^ 



(87) 



Parametrizing the Fermi parabola in the usual way, z(s) = {^/vn + is) 2 , and using 
(|87p we obtain 



■ tr Apf(z) dz 



< 



— tr An(z(s)) ds 



4 yW+T 



_, =ds- 

m. (vn + s 2 )2 ^/\s\ Jr {v N + S 2 )yJ \s\ 



Co 



ds 



Because the integral 



( 7= ds = 4A 

Jn(v N + s 2 Y^ 



o (l + s 4 ) a 



ds 



exists in particular for a £ {fjl} we conclude that the integral in (|55|) converges 



absolutely and satisfies the bound given there with an appropriate constant. 



□ 



5.2. Dominant term 



To begin with, we single out the logarithmic part of the integral over tr Bn(z). 

Proposition 5.2. Let V G £ 1 (R) satisfy (I5T1) . Furthermore, assume that V satisfies 
(1291) and (|30p with some a > and < p < 1. The following integral converges 
absolutely and behaves asymptotically for large N and L with p = (N + |)/(2i) as 

1 f 1 

— / -tri?Ar(z)dz = xjv7l +£jv (88) 
Z7rl Jr N z 

with x/v from Provosition \3.6\ and the coefficient "/l := j s .l + 7c. L, 

j s , c , L := - 1 _ / ^-(Sl(z( S )) Ws!C (^),f(^)!l(z( S )) Us , c (^))&, (89) 
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where z(s) — (^/vn + is) 2 - The operator 

F{v N ) := J* y/\V\(P s (u N ) + P c (v N ))y/\V\J. 
The error term can be estimated as 

\e N \ = 0(1 + — lniV + — In 2 N). 



(90) 



(91) 



Proof. We proceed in three steps. We show that the integral converges absolutely and 
then apply the delta-estimate, Proposition 14. 11 as often as possible. Finally, we single 
out the dominant part of the truncated resolvent. 

(a) To begin with, we bound tr P>n(z) for z € Tjv, i.e. yfz = ^vn + is. With the aid 
of Lemmas 13.11 13.21 and 13.51 we infer from Formula (1841) 



(92) 



ixB N (z)\ < \\n(z)\\\\n(z)\\\\y/\V\S N (z) V \V\\\x 

X (\d s (z)\\\uj s (z)\\\\Us(z)\\ + |do(2!)|||w c (z)||||w c (2!)|| 

<*c D clcl\\v\w^ 1 

IT 



f w L ,*( a ) ln(iV + l). 



Parametrizing the Fermi parabola as usual, z(s) — {y/vji + is) 2 , we conclude that 

1 



1 f 1 1 
— / -tr B N (z)dz = - 



2ni 



7T J u y/VN + IS 



tr P>n(z(s)) ds 



(93) 



converges absolutely because of Lemma 14.21 
(b) For the sake of brevity, we define 

&s,l(s) ■= d s (z(s)), S C:L (s) := d c (z(s)). 



(94) 



The following calculations look alike for 5 S! l(s) and S c _l{s) and the corresponding 
quantities. Therefore, we simply write 8l(s) etc. to denote either case. We evaluate 
the integral in (|93[) and start with replacing Sn(z) by Sn{^n)- This causes an error 



B. \/ v N + IS 

x (Jfl(z(s))uj(z(s)), y/\V\(S N (z(s)) - S N (v N ))^\V\jn(z{s))uj(z(s))) ds 



< c D clcl 



L 



Recalling Lemma T3. 5 1 we use Proposition 14. II with 
h(s) = \\V\V\(Sn(z(s)) - S N (u N )) 



w L , a (3)\\y/\V\(S N (z( 3 )) - Sjf{y N )W\V\\\ ds 



§{L) = InL, Q(L) = L 
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and obtain the error 



= (1) <C 2 fl 



We replace w(z) by w(i/jv). This leads to 

|(jo(^( s )) w (z(s)),yjy|5 JV K)v^ 7 l^(^(s))^(2(s))) 
- (m(z(s))u(p N ), ^s N (p N )^\m(z(s))uj(u N ))\ 

= | Re(m(z(s))(u(z(s)) - uM), y/\V\S N (u N )y/\V\JQ(z{a))(u(z(a)) + u(v N )))\ 
< \\n(z(s))\\^u(z(s)) - u( m )\\(\\uj(z( S ))\\ + \MvN)\\)\\\\y/\V\S N (v N )y/\V\\\. 

To bound the error by dint of Proposition l4~Tl we recall Lemma l3~2l Actually, one could 
estimate ||cl>(^v)|| via (|46D . However, it also possible to use the bound for u>(z(s)) which 
simplifies the calculations. Then, 

f L (s) = e- 2Ls \\cj(z(s)) - u(y N )\\\\u(z{s)) + u(y N )l 0(L) = 1, 0(i) = L 1 ^, 

where the rate functions are from Lemma 13.21 Thus we get 



e^<C 3 (^ + (l-p) l ^)lnN, 



where the rightmost logarithm is from Lemma 13.51 

(c) We decompose Sn(vn) according to Proposition 13.61 and hnd 

\\VW\S 3 ,n(vn)W\\\ <C, j = 1,2, \\V\V\Si,n(vn)VW\\\ < C 
because of Lemma [A. 2 1 Thus, we are left with the integral 

5l{s) ^{z{s))u,{z{s))fds<C D CiCtL [ WL ^ s \ ds<C(l + ^) 



E 



^/v N + s 2 7r Vvn + s 



where we used Lemma [4.2l Hence, the remaining term is given via Sq,n(vn)- Rewriting 
it by dint of the addition theorem for the cosine yields the statement. □ 



To determine the asymptotics we have to study the coefficients 7^. 

Theorem 5.3. Let V € X^R) n L°°(R) be non-positive and X 2 V € L 1 (R). Assume 
further that V satisfies (|61[) . (|64p . (|69l) . and the conditions (|29p . (|30p with some < 
p < 1 and a > 0. In £/ie thermodynamic limit according to (|80p . </ie Anderson integral 
has the leading asymptotics 

1 N = 7 In TV + o(ln TV) (95) 

wif/i £/ie constant 

7:=8^to[(l + ^, )-^ >( J>0 (96) 
and £/ie 2x2 matrix &oo,Q as i n (|68p . Moreover, if ||V||i is small enough then 7 > 0. 
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5. Asymptotics 



Proof, (a) We use (f89|) and show that jl 7- To this end, we substitute s = t/L in 
(l89l) and obtain for "f s> L the integral 



TR. 



(V^v + f) shr (L^7 + it) 



x (n((v^ - |) 2 K(^), F(u N )Sl((y/UH + *j-) 2 )u.(v N )) dt, (97) 

where the minus sign conies from the complex conjugate. We want to use Lebesgue's 
dominated convergence theorem. To begin with, the integrand has the integrable 
bound 



1 



-.C D Cl\\F{y N )\\ 



by Lemmas 13.11 and 13.81 The convergence of the first factor in (1UT)) is obvious. From 
Lemma 13.91 we know 

||fi± {u N , ±t) - n((V^± |) 2 )|| < C' n J-^(l + L^wl^j)^) 

for all t € H. Here the plus sign is used for even N in z/jv and the minus sign for odd 
N . For each subsequence we have convergence for fixed t. We first look at the limit 
7s, l — > for even N. Here, 



1 



7 



1 



87r 2 ^j\r J sr. (coshi + i sinhi) 2 
where we used (f3"9"|) for the sine. The same reasoning yields 



(i/jy, -t)uj s (v N ), F^n)^ (v N ,i)u; s {v N ))dt, 



7, 



r (n+ (^Ar,-t)w c (l/jv),-F 1 (^Jv)^i,(^A',O a; c(^Jv)) dt. 



8tt 2 vn (cosht — i sinht) 2 °° 
(b) To evaluate 7 we recall the definition (j90|) of F(un) and obtain 
(fi+ (i/jv, -t)w a , c (i4v),.F(i/jv-) fi £ (^iv,i)w S)C (^jv)) 

+ (Jfi+ (i/jv, -t)w a , c (i/ w ), y/W\Pc(m)VWJtt£(j'N,t)u l ,, c (v N )) 

= (Jft + (v N , -t)u atC (v N ),U a (l/ N ))(u a {l/ N )JQ,^(uN,t)cj a>c {VN)) 

+ (y N , -t)u} StC {vN),0J D {v N )){w c {v N )JQ^{v N ,t)ui 8)C {v N )). 

Therefore, by using the matrices f2+,(i/jv,t) and f(t) (cf. (|67D and (|71|l) we can write 
the integrand of 7+ + 7+ as the trace of 2 x 2 matrices which leads to the integral 

I-=~h I trfr'^fi+K,-*)*^^)] dt 



1 

2/ 



TR. 



tr [T l {t)Vt oofi {i' N )Z + {v N , -t)*Z + (^ A r,t)f2 OOi0 (^jv)] dt. 
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Here we used (1721 to express f2+ (fN,t) through ^oo.o^aO and Z + (vN,t). With the 
aid of the relation 

sXo,o(^)^ + (^,^)* = Z+^n^^o^n) (98) 
and the cyclicity of the trace we obtain 

1 = -— I ti[Z + (v N ,t)h oofi (v N )f\t)Z + (v N ,t)Vl oofl (v N )\dt 



2i 



R 



tr [(2' + )'(^Ar,t)r2 00! o(^Ar)] 



= Um tr [(Z+(v N ,t) - Z+{v N , -t))Cl^ ]. 

2,1 t-foo J 

We compute the difference 

Z+(u N ,t) - Z+(y N ,-t) = 1 r^=Z+(v N ,t)n oo , (vN)(T{t) - T(-t))Z+(v N ,-t) 

2iImT(t) A , . , - , , A , . 

V 7 -Z+(^ A r,t)r! OOi0 ( i y A r)Z + (^ A r,-t). 



Thus, our integral becomes 

I= - lim Im{r(i)tr [Z + (^at, t)^^ (^Ar)Z + (^Ar, -t)^^ (^at)1 } 
- - lim Imr(t)tr [Z + (i/jv, t)^oo.o(^) 2 ^ + (^, 0*] , 
where we used again. Since we have (see (|4TJ1) and ([71])) 

lim Rer(i) = 0, lim Imr(t) = — 1, lim f(i) = —it, 

t—±00 t— fOO t— S-OO 



we obtain 



/ = tr [(1 - — — noo,o(fjv)) _1 (l + TT^^oo.O^Jv)) ^oo.Q^iv) 2 ] 



= tr [(1 + -^-^oo.o(^A') 2 ) _1 ^oo.o(^) 2 ] • 

Apart from the prefactor this is the cocfhcicnt 7 for even TV. For odd N analogous 
computations lead to the same 7. 

(c) From r2oo,o(K/v)* = ^oo.o^at) we infer that 7 is the trace of the product of two non- 
negative matrices. Hence 7 > 0. Toshow7 > first note that H^oo.o^aOII = 0( || V|| 1 ) . 
Thus, 

7 = tr(lW^) 2 )+0(||^||i). 
Furthermore, with the aid of Lemma 13.91 we obtain 



linoo.oCi/jv) - 1|| < ^\\toaoA v *)\\\\yfiv\ G My/W\\\ < -^c nao0 \\v\\ u 
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which implies for the matrix elements 

(W s , c (z/jv),^oo,o(^iv)w s ,c(^jv)) = {0J S: c{vn),^s,c( v n)) + 0{\\V\\l). 

We conclude that 

j = tr(fll) + 0(\\V\\l) 

with the matrix 

^ {{u s {vn),u s {v n )) {<jJ s {vn),u c {v n ))\ 
1 ' \(wc(vn),u s (v n )) (uj c (v n ),u c (v n ))) ' 

For simplicity let us estimate 

tr(0 2 ) > -(tr^) 2 = \{\\lo s { Vn )\\ 2 + |KK)|| 2 ) 2 - 
Writing out the scalar products we obtain 

IMz'aOII 2 + ||w s Oat)|| 2 = / \V(x)\(sin 2 (u N x) + cos 2 (i/jva;)) da; = \\V\\ 1} 

Jul 

which finally shows 

1>\\\V\\1 + 0(\\V\\1)- 
Thus, 7 > for ||V||i small enough. 

(d) The integral formula (15TT) involves Hm and therefore differs from the actual An- 
derson integral by 

\tiP N (i-n N )-trP N (t-n M )\ = |trP w (n w -n M )| < |jv-m|. 

Using V(x) < and Proposition 13 . 1 II we conclude 

N<M<N+- + C E 

with E — vjs[. Therefore, replacing II by II m causes an error that can be bounded 
by a constant. □ 

Some remarks seem to be appropriate. 

1. The potential being non-positive is nowhere needed except for the final step in 
the proof of Theorem 15.31 It is only to ensure that the integral formula (|5T|) 
indeed represents the Anderson integral correctly. 

2. Potentials with compact support do not require any sign conditions. In that case 
one can avoid the Bargmann inequality in Proposition 13.111 and prove estimates 
for the eigenvalues by using square well potentials instead. 
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3. One could have used Lebesgue's theorem already in Proposition 15.21 instead of 
the delta-estimate |4~T1 and thereby treated the integral in one step. However, the 
delta-estimate gives explicit error bounds instead of the mere o(ln N) from the 
dominated convergence. 

4. The coefficient 7 in can be cast into another form 



Finally, we obtain an upper bound on the transition probability, T>n — \(<f> , ip N )\ 2 , 
between the two ground states <fi and ip . 



Corollary 5.4. Under the conditions of Theorem \5.3\ and assuming that the coefficient 
7 is strictly positive one has for any e > and large N the asymptotic bound 

T> N < e- lN < N-~< +e 

on the transition probability T>n (see (HJj. 

Proof. We follow Anderson's arguments and use Hadamard's and Bessel's inequality 
to obtain 

N N 

|det((^,^))| 2 <nEl(^'^)| 2 
j=l k=l 

N 00 

j = l k=N+l 
N 00 

j=l k=N+l 

= exp(-Zjv) 

where we applied ln(l — x) < —x in the last step. If we choose N large enough then 
the error satisfies (see (|95)l ) 

gy-7lnJV| 



In TV 

Hence, for large N 



< e 



V N < cxp(-Ijv) < exp(-7lniV+ 7 " — —hxN) < exp(-flnN + elniV), 

In AT 

which proves the statement. □ 
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A. Estimates 



At various points we need estimates which are not directly related to our main subject. 
To begin with, we discuss the following sums 

N N-l N N 

V - = 2 Q V - < 4 a V - < 4" / - dt 

j^{N + \-j) a ^( 2 J' + 1 ) a " + Jo (t + l) a 

To evaluate the integral we distinguish the cases < a < 1, a = 1, and a > 1. We 
obtain 

n f^rr for a > 1, 

E7^-T— v^<4 Q T^^ + l) 1 - tea<l, (99) 
J=i 1 2 7j [ln(7V + l) fora = l. 

The constant in Proposition 13.61 requires more reasoning. 
Lemma A.l. Let M £ R, w £ (D smc/i £/ta£ Rew + i - M > 0. Then, 



e 



wt cosh(Mt) ^ < 2 

cosh | ~~ «; + | - |M|' 



In particular, \Zcn\ < 4 m Provosition \3.b\ Furthermore 
_ t sinh(Mt) 



sinh | 



■fif 



M , w + | + M 

1 ln- 



(w + i) 2 — M 2 w + \ - M 

1 f°° y 1 

+ 8M(w + -)J^ — - | - — - — — - j _ — - — ^ _ - dy. 

In particular, this yields the asymptotics for all N G N 



3G 



M .sinh((A+ ±)<) , N 2A + 1 
-(iv+ 7 )t " ut rft = ln(4A + 3) + ^y + CAr withO<c N < 



2N+1 

e 



o sinh | v ' 2A^ + § - 2N + f 



Proof, (a) The first estimate follows from estimating the integrand via 

-e* < coshf < e*, t> 0. 
2 - - ' - 



(b) To begin with, we write for t > 



1 i 00 



sinh | 1 — e 

2 j=0 
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dt = 2M ' 



1 



and integrate termwise which yields 
f°° c _ wt smh(Mt) 

Jo sinh| ""^(«,+ i+j)a_Afa- 

We apply the Abel-Plana summation formula [6, Th. 4.9c] 

2M 



sinh | 



1 



1 

2 {w + ±) 2 - Af 2 ./„ 
1 



1 



(to + | + x) 2 - M 2 

1 



dx 



{w + \+ iy) 2 — M 2 (to + \ - iy) 2 - M 2 



3 2try _ I 



dy 



which implies the formula. Finally, 

< 

yields the estimate. 



< 



;27ry _ 1 - 27T 



□ 



It is obvious that the kernels So ; pf(z;x,y) and So > N(z;x,y) in Proposition 13.61 are 
bounded functions in x and y as long as s/z € 1R. For the other two, S\,n(z; x, y) and 
Sx t ]\f(z; x, y), this is only true at the Fermi energy and requires some computation. 

Lemma A. 2. There is a constant < C < oo such that for all a € 1R and N 6 N 



• //„ In/ nnC°s((^+ §)«) , 

sm((JV + -)(«- a)) » d« 

q Z COS 2" 



< c. 



Proof. By a change of variables we may assume a > 0. For brevity, we put M := iV + 1 
and a := [^"J ■ Using the periodicity of the integrand we get 

sin(M(u-o)) — —n-^du = (1 - (-1 Q / sin M tt-a —^du 

sm 2 ^0 Sm 2 



+ (-lf 



To begin with, let < b < vr. Then, 



sin(Af (u — a)) 



sin(Mu) 



sin(M (tt — a)) 



sin(Afu) 



<iii. 



du = 2 sin(M(w - a)) 



sin u 



sin(M(tt — a)) 



-sin(Afti) sin(Afu)- 
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A. Estimates 



Via the Taylor expansion of the sine one can show that the rightmost integral is 
bounded in M and a. A change of variables in the remaining integral yields 

b ■ tin ..sin(Mu) f Mb sin it 

sm(M(u — a)) du = I smii-Afo du, 

u J u 

which can be shown to be bounded in M and a via integration by parts. The case 
7T < b < 2tt can be reduced to this one. The cosine integral can be treated in like 
manner. □ 



Lemma A. 3. Let WL, a be as in (|48|) with a > 0. Then for all s > and L > 



Proof. We divide the domain of integration and note a > 0, 

L -, r L/2 i r L 

-e 



For Ls > 1 we have 



whereas for Ls < 1, 



2sx dx= / -t— e 2sx dx + / -—e 2sx dx 



roo i ol+« i 

<e Ls - dx + - -(e 2sL -e sL ) 

Se J (i +a: )i+a ax+ L ^2s ( ' 



e 2Ls (-e- Ls + -^(l-e- Ls ) 



(\- e - Ls \ < -=-1 



2 a ^ _ 2 a 2 a s a 



/ 1 — i/S\ ^ 

L!+«s l j - L a ~ (Ls)* 1 

where we used 1 — e~ x < x for x > 0. This concludes the proof. □ 
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